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Abstract 
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■ When the strong CP problem is solved by spontaneous breaking of an anoma- 

\& , lous global symmetry in theories with gauge-mediated supersymmetry break- 

| ing, the pseudo Goldstone fermion (the axino) is a good candidate of a light 

J> ' sterile neutrino. Its mixing with neutrinos relevant for current neutrino ex- 

periments can arise in the presence of R-parity violation. The realistic four 
. neutrino mass matrix is obtained when the see-saw mechanism is brought in, 

and an ansatz for the right-handed neutrino mass is constructed. 
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Current neutrino experiments observing the deficits in the atmospheric [|IJ and solar 
neutrino fluxes possibly hint at the existence of an SU(3) C x SU(2)l x U(1)y singlet fermion 
that mixes with the ordinary neutrinos. More excitingly, the reconciliation of the above 
neutrino data with the candidate events for the — * v e oscillation || requires the existence 
of such a singlet fermion (called a sterile neutrino [[|) with the mass ^ (9(1) eV. Since the 
mass of a singlet state cannot be protected by the gauge symmetry of the standard model, 
the introduction of a sterile neutrino must come with a theoretical justification. In view 
of this situation, there have been many attempts to seek for the origin of a light singlet 
fermion with the required properties [|5],|Bj . In this letter, we point out that a sterile neutrino 
can arise naturally in a well-motivated extension of the standard model, namely, in the 
super symmetric standard model (SSM) incorporating the Peccei-Quinn (PQ) mechanism 
for the resolution of the strong CP problem 0, and the mechanism of supersymmetry 
breaking through gauge mediation ||. 

Most attractive solution to the strong CP problem would be the PQ mechanism which 
introduces a QCD-anomalous global symmetry (called the PQ symmetry) spontaneously 
broken at a high scale f a ~ 10 10 — 10 12 GeV, and thus predicts a pseudo Goldstone boson 
(the axion, a) [[?[]. In supersymmetric theories, the fermionic partner of the axion (the axino, 
a) exists and would be massless if supersymmetry is conserved. But supersymmetry is broken 
in reality and there would be a large mass splitting between the axion and the axino, which 
depends on the mechanism of supersymmetry breaking. In the context of supergravity where 
supersymmetry breaking is mediated at the Planck scale Mp, the axino mass is generically of 
the order of the gravitino mass, 1713/2 ~ 10 2 GeV 0, which characterizes the supersymmetry 
breaking scale of the SSM sector. However, the axino can be very light if supersymmetry 
breaking occurs below the scale of the PQ symmetry breaking, as in theories with the gauge- 
mediated supersymmetry breaking (GMSB) ||. The GMSB models are usually composed of 
three sectors: the SSM, messenger and hidden sector. The messenger sector contains extra 
vector-like quarks and leptons which conveys supersymmetry breaking from the hidden sector 
to the SSM sector 0. 



In order to estimate the axino mass in the framework of GMSB, it is useful to invoke the 
non-linearly realized Lagrangian for the axion superfield. Below the PQ symmetry breaking 
scale f a , the couplings of the axion superfield $ can be rotated away from the superpotential 
and are encoded in the Kahler potential as follows: 

I I Ja 

+ higher order terms in /„ (1) 

where Cj is a superfield in the SSM, messenger, or hidden sector, and xj is its PQ charge. 
Upon supersymmetry breaking, the axino gets the mass from the third term in Eq. ([[]), 

mai ~ Xjy- (2) 

J a 

where Fj is the F-term of the field Cj. If there is a massless fermion charged under the 
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PQ symmetry, the axino would have a Dirac mass of order Fj/f a ||1 1||P|. It is however 
expected that there is no massless mode [except neutrinos, see below] in the theory, and 
each component of the superfield Cj has the mass of order Mj ~ y/Fj unless one introduces 
extra symmetries to ensure the existence of massless modes. Then, the axino mass is see-saw 
reduced to have the majorana mass: 



vn-n 



( 3 ) 



Now one can think of three possible scenarios implementing the PQ symmetry: PQ symme- 
try acting (i) only on the SSM sector, (ii) on the messenger sector, (iii) on the hidden sector. 
In each case, y/FJ is of order of (i) 10 2 - 10 3 GeV, (ii) 10 4 - 10 5 GeV, (iii) 10 5 - 10 6 GeV. 
In the last case, we restricted ourselves to have a light gravitino m 3 / 2 < 1 keV which evades 
overclosure of the universe if no entropy dumping occurs after supersymmetry breaking. 
Then we find the following ranges of the axino mass: 

(10- 9 -10- 6 )eV(i^^) 2 (i) 

(10- 3 -l)eV(±2^) 2 (ii) (4) 



;i_l 3)eV(i^^J (iii) 

Interestingly, the axino mass in the case (i) or (ii) falls into the region relevant for the current 
neutrino experiments. Our next question is then how the axino-neutrino mixing arises. 

The mixing of the axino with neutrinos can come from the axion coupling to the lepton 
doublet L, x u (& + $)L^L, which yields the mixing mass, 

F v 

m au ~ X u — (5) 
la 

where x v is the PQ charge of the lepton doublet. It is crucial for us to observe that nonzero 
F v can arise when R-parity is not imposed. To calculate the size of F u , we work in the basis 
where the Kahler potential takes the canonical form. That is, L'H\ + h.c. is rotated away 
in the Kahler potential, and the superpotential of the SSM allows for the bilinear terms, 

W = + e l fiL i H 2 . (6) 

where the dimensionless parameter measures the amount of R-parity violation. Due to 
the e-term in Eq. (|6|), one has F v = e^v sin/? where v = 174 GeV and tan/3 = (H 2 )/(Hi), 
and therefore, 

4 ( 6 t sin/3 \ /KPGeV\ ( p \ 
m -~ 10 eV (27l^)(^^)(v300GevJ- (?) 



1 Note that the late-decaying particle scenario for the structure formation with the axino in the 
MeV region and the gravitino in the eV region [12] can be realized in this case for F ~ (10 5 GeV) 2 . 
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Since R-parity and lepton number violating bilinear operators (0) are introduced, we have to 



take into account the so-called tree- level neutrino mass |13| . The tree mass arises due to the 
misalignment of the sneutrino vaccum expectation values with the £j terms, which vanishes 
at the mediation scale M m of supersymmetry breaking, but is generated at the weak scale 
through renomalization group (RG) evolution of supersymmetry breaking parameters. This 
tree mass takes the form of m}J ee oc e^ej, and its size (dominated by one component ei) is 
given by E , I5 |, 



m* ree w 1 eV 



3 x 10-V \My 



Shi, M„ 



where a f ~ e f sin£ { — T ) { — In — ■ J (8) 



where Mz, Mi/ 2 and mj are the Z-boson, gaugino and slepton mass, respectively, and hj,, At, 
are the 6-quark Yukawa coupling and corresponding trilinear soft-parameter, respectively 
In Eq. (g), the term proportional to /i 2 In (M m /m z ~) characterizes the size of the RG- induced 
misalignment. Taking tan/3 = l,/i.A& = m~, and M m = 10 3 m z ~ as reference values, one finds 



m 



,tree i n— 4 



\io-y \m 1/2/ 

which grows roughly as tan 4 (3 in large tan (3 region. 

Based on Eqs. (f|), (0) and (|9|), we can make the following observations: 

• The just-so solution of the solar neutrino problem with large mixing and Am 2 ol m 
0.7 x lCT 10 eV 2 implies m* r e ee ,m s < ma„ 8 ~ 10~ 5 eV, and could be realized in the 
case (i). For this, we need f a ^ 10 11 GeV and e\ ~ 1CT 7 . 

• The small mixing MSW solution requiring # so i «4x 1CT 4 and Am^ ps 5 x 1CT 6 eV 2 |2| 
can be realized for the case (i) or (ii) if < rria ~ 2 x 10~ 3 eV and ps 10~ 4 eV. 
For this, one needs f a ~ 10 10 GeV and tj ~ 10~ 8 in the case (i); or f a ps 10 12 GeV and 
£j ~ 10~ 6 in the case (ii). 

• The — > a explanation of the atmospheric neutrino oscillation |16| requiring nearly 
maximal mixing || is realized if rria Ue > ma,vn^™ e and Am 2 tm ps 2ma Ve { m a + m !^ ee ) ~ 
3 x 10~ 3 eV 2 . The best region of parameters for this is e 2 ~ 10~ 5 and f a ~ 10 10 GeV 
prefering the case (i). 

Note that a low tan (3 is preferred to suppress m* r . ee in all of the above cases. 

Having seen that the v e ^ — a mixing can arise to explain the solar or atmospheric neutrino 
problem, let us now discuss how all the experimental data [|l|-|3| can be accommodated in 
our scheme. Recall that there exist only two patterns of neutrino mass-squared differences 
compatible with the results of all the experiments. Namely, four neutrino masses are divided 



into two pairs of almost degenerate masses separated by a gap of y Am1 SND ~ 1 eV as 
indicated by the result of the LSND experiments [[3], and follow either of the following 
patterns |I7J]: 
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atm solar 

(A) mi < m 2 <C m 3 < m 4 , 

V v ' 

, LSND 
solar atm 

(B) 777,1 < 777.2 <C 7773 < 7774 . 

LSND 

In (A), Am^ is relevant for the explanation of the atmospheric neutrino anomaly and 
Ar/743 is relevant for the suppression of solar i/ e 's. In (B), the roles of Am 21 and Am^ are 
interchanged. 

In our scheme, the required degeneracy of m 3 and m 4 could be a consequence of the quasi 
Dirac structure: m aVi 3> m^, 777* r . ee . To have 7774,7773 « ~ leV, we need a large value 
of €i ~ 10~ 4 e.g. for f a = 10 10 GeV. But this makes the tree mass too large (777* ree ~ 1 eV) 
to accommodate A777 2 ol or Am 2 tm for (A) or (B), respectively. For the pattern (B), the 
atmospheric neutrino oscillation could also be explained by the — v T degeneracy with the 
splitting 7774 — 7773 ~ 10~ 3 eV. However, when the ordinary neutrino mass comes from R- 
parity violation, the neutrino mass takes generically the hierarchical structure since the tree 
mass (of the form m^ 66 oc €i€j) is rank-one. Even though this structure can be changed due 



to the one-loop contribution through the squark and slepton exchanges p^3|JT5| , one needs 
fine-tuning to achieve the required degeneracy: (7774 — 777 3 )/777 3 ~ 1CT 3 fT5|| . 

The simplest way to get the realistic neutrino mass matrix would be to introduce heavy 
right-handed neutrinos iV whose masses are related naturally to the PQ scale f a ~ 10 12 GeV. 
For this purpose, let us assign the following U{1) PQ charges to the fields: 



(10) 



H 1 H 2 L N (j) (j)' a a' Y 
1 12-3-116-60 

Then the PQ symmetry allows for the superpotential W = Wn + Wpq where 

M P Mf, (H 2 ) 2(a) 

W PQ = (<M>' + aa' + Ml)Y + M 2 Y 2 + V 3 + ~«/>V + ~0'V (11) 

3 6 6 

where My, Mi, M 2 ~ f a . Minimization of the scalar potential coming from Wpq leads to 
the supersymmetric minimum with the vaccum expectation values (0), (</>'), (a), (a') ~ f a 
and (Y) ~ satisfying the conditions (0/0') 6 = (a'/a) and (4>4>') = 6 (era'). Note that the 
superpotential Wn provides a solution to the /i problem Jig] , and generates the right value 
of €i for our purpose, that is, 

V|^10 2 GeV, ei ~-^^10- 6 . (12) 
P M P ' h^Mp v ' 

Now the usual arbitrariness comes in the determination of the Dirac neutrino mass mf and 
the right-handed neutrino mass M^. Guided by the unification spirit, let us first take the 
Dirac masses, mf ~ 1 MeV, mf ~ 1 GeV and mf ~ 100 GeV following the hierarchical 
structure of up-type quarks. Then, the see-saw suppressed neutrino mass matrix 7?7 3x3 = 
—mfmfM^ 1 takes the form, 
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m 3x3 ~ 1 eV 



/ 10- 8 c n 10^ 5 ci 2 10- 3 ci 3 
10- 5 c 12 10- 2 c 22 c 23 i (13) 



V 10- 3 c 13 c 23 10 2 c 33 

where = M^/M^ 1 and m^mfM^s ~ leV. Note that the matrix ( [IB] ) is close to the 
neutrino mass matrix ansatz ||19|| : 

/ 0.01 \ 

(14) 





( o 





m 3x3 ~ 1 eV 





io- 




V o.oi 


1 




which explains the atmospheric neutrino and LSND data, except that (|l3l) has a too large 
component, wi 33 3 . Restricted by minimality condition and requiring c 33 = 0, the ansatz ( |14"D 
can be translated to the corresponding one for the right-handed neutrino mass matrix: 

(A B C\ 
M i:j = B B 2 /A . (15) 
\C J 

The parameters A, B, C in Eq. (|TJ|) are determined from the observational quantities as 
follows: 



— — Am^ m 3 



^LSND - (16) 



AmL _ C m$ 



Am LSND B m 3 

from which one finds A, C ~ 0.15, and B 2 / A ~ 105 for B w 10 11 GeV. 

In conclusion, it has been shown that the axino, which is predicted by the PQ mechanism 
in supersymmetric theories, can be naturally light and mix with ordinary neutrinos to explain 
the solar or atmospheric neutrino problem in the context of gauge mediated supersymmetry 
breaking. The lightness is ensured by the small supersymmetry breaking scale y~F ^ 10 5 
GeV (far below the PQ symmetry breaking scale f a ) and the mixing is induced due to the 
presence of the R-parity violating bilinear term ej/iLjif 2 . The v e — a mixing can arise to 
explain the solar neutrino deficit when the parameters are in the ranges: e% ~ 10~ 6 and 
f a ~ 10 12 GeV for which the axion can provide cold dark matter of the universe. To account 
for the atmospheric neutrino oscillation [y^ — > a), larger e and smaller f a are preferred: 
e 2 ~ 10" 5 ,/ a ~ 10 10 GeV. 

With R-parity violation alone, it is hard to find the four neutrino mass matrix which 
accommodates all the neutrino data. It has been therefore attempted to obtain a realistic 
four neutrino oscillation scheme relying on the see-saw mechanism combined with R-parity 
violation. In this case, only the pattern (B) with the u e — > a solar and — > v T atmospheric 
oscillations can be realized. A toy model has been built in a way that the right values of \i 
and e are generated naturally through the PQ symmetry selection rule. In the unification 
scheme where the Dirac neutrino mass follows the hierarchical structure of up-type quarks, 
an ansatz for the right-handed neutrino mass matrix has been constructed to reproduce the 
required light neutrino mass matrix. 
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